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Bulk viscosity has been intrinsically existing in the observational cosmos evolution with various 
effects for different cosmological evolution stages endowed with complicated cosmic media. Normally 
in the idealized "standard cosmology" the physical viscosity effect is often negligent in some extent 
by assumptions, except for galaxies formation and evolution or the like astro-physics phenomena. 
Actually we have not fully understood the physical origin and effects of cosmic viscosity, including its 
practical functions for the universe evolution in reality. In this present article we extend the concept 
of temperature related viscosity from classical statistical physics to observational cosmology, espe- 
cially we examine the cosmological effects with possibility of the existence for two kinds of viscosity 
forms, which are described by the Chapman's relation and Sutherland's formula respectively. With 
that the unified dark fluid model and also a modification of standard model with viscosity which 
is named as ACDM-V model are constructed physically. While the unified dark fluid model con- 
strained by observational data is closely related with bouncing phenomenon, the ACDM-V model is 
supported more by data-sets fittings with results indicating that temperature related viscosity could 
not totally explain the rich physics of dark energy phenomena. In addition to the enhancement to 
cosmic age value, the ACDM-V model possesses another two pleasing features: the prediction about 
the no-rip/singularity future and the mechanism of smooth transition from imperfect cosmological 
models to perfect ones. 



I. INTRODUCTION 

The perfect fluid cosmological models are built on the dynamics governed by Einstein's theory of general relativity 
with the cosmological principle validity at large cosmos scales and the assumption of cosmic compositions modeling as 
an idealized perfect fluid, which means that all components of the matter-energy in our universe are considered as per- 
fect fluid without any viscosity or which may be interpreted as a result of having taken the gross grain approximation. 
The well-known ACDM model, which is regarded as the most idealized standard description for current cosmology, 
obtains a constant dark energy term (the famous cosmology constant) that permeates the space everywhere and glob- 
ally fits the observational data sets very well (2{|. But researches have shown that constant dark energy models are 
not well confirmed by both observations and theoretical considerations [H, • So far the related research work in this 
line has produced a variety of improved models, and one direction is studying the practical effects of the physically 
existing viscosity on the cosmic evolutions and its manifestations 0, [l2| • It is a natural fact that viscosity exists in 
the universe evolution history with its eminent role in the thermal plasma as from hot Big Bang like stages and the 
important effect of turbulence viscosity has been observed both in intra-galactic and inter-cluster medium. 

Even though the existence of viscosity is not a problem, what has confused us is its possible origin in the cosmos 
and how the cosmic viscosity affects the cosmological evolution. In other words, we are searching after physically 
reasonable explanations and acceptable behaviors of the large-scale existing viscosity in the universe. 

It is easy to break down the assumption of perfect fluid by proper modification, but not the cosmological principle 
which is not only a theoretical principle for large scale cosmology but also an observational feature of our cosmic 
microwave radiation background. There are two methods for conceiving imperfect fluid models: one obeys the 
cosmological principle which means the momentum-energy tensor is diagonal in the diagonal metric compatible with 
the cosmical principle; the other violates the cosmological principle in local region but should be subjected to an 
average procedure, so that after which the statistical homogeneous and isotropic properties are retrieved at large 
scales. The first method is widely followed in previous literature, while the second one which is also called back- 
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reaction is famous for its complication and is still in progresses Q. 

In this paper, the first method will be adopted in order to build temperature related imperfect cosmological models 
rather directly. So we arrange the contents as: In the following sub section, before introducing the model buildings 
we simply review some basic concepts which are closely related with our theme in this present work. In the second 
section, we will discuss two kinds of imperfect fluid models with a temperature related viscosity origin. The first 
kind is an unified model where matter and dark energy are indistinguishable clearly or coupled each other, while 
the second is a modification of the ACDM model where the bulk viscosity acts as an extra pressure of total fluid in 
addition to the normal one. Model constraints with cosmic observational data sets fittings are also included in the 
second section. We discuss and perform some numerical contrasts for the ACDM-V model to ACDM model in section 
three. In the end of this present work, we briefly discuss some features of the new models numerically with the small 
deviations away from the ACDM model and give some conclusions in the fourth section accordingly. For the benefit 
of the related community we present in details how to do the statistic analysis for the cosmology model buildings with 
the currently available observational data-sets from direct measuring H(z), SNe la from the Union 2.1 and baryon 
acoustic oscillation (BAO) from large scale surveying as the SDSS (BOSS data-sets) and WMAP-9yr (including the 
6dFGS [H] (Beutler et al. 2011), Wiggle Z [H] (Blake et aJ.2011), and SDSS III survey results) in the appendix. 

A. Einstein's equations 

Every homogeneous and isotropic cosmology model is born in the environment of Einstein's dynamic governing 
equations in the Friedmann Robertson Walker (FRW) metric [HIHII, which read (Since the curvature of the universe 
is theoretically predicted and observed today to be rather close to zero, so we focus on the flat metric geometry in 
this article). 

H 2 = n 2 p (1) 

H 2 +H=- = -^-(p + 3p) (2) 
a 2 

where H = - represents the Hubble parameter, alternatively, we usually use E = H/H (the dimensionless Hubble 
parameter) for convenience. Notice that, H is the observed Hubble parameter at present. The combination of the 
Einstein's component equations (Friedmann equations above) gives the continuous equation 

p + 3H(p + p) =0 (3) 

which means as long as we follow the Einstein's equations, the total energy is automatically conserved. Notice that 
effective (or total) density p and pressure p are concerned in Eqs. (H|)— (J3J) . 

In addition to the basic functions above, we also pay attention to the expression of equation of state (EOS) parameter 
which is of vital importance in distinguishing phases of evolution of the universe: 

w = - = 7 (p) - 1 (4) 
P 

In the well-known ACDM model, w can be reduced to —Uj^/E 2 according to the assumption that radiation occupies 
so little fraction of the total density that could be neglected in the late stage evolution. 

B. viscosity and its classification 

Generally speaking, viscosity is still a phenomenological concept generated from the inhomogeneous motion and 
molecular- like interaction in the fluid [33||. Mathematically the classical viscous stress tensor could be expressed as 

cr V is;ij = C,d k v k 5 i: j + r](diVj + djut - ^d k v k Sij) (5) 

where v k is the k-dimcnsional component of velocity. 

According to the theoretically homogeneous and isotropic property of the large scale observational universe discussed 
above, we focus on the first term in Eq. namely the bulk viscosity (or the 2nd viscosity) term, while neglecting 
dissipation and shear viscosity which are incompatible with the cosmological principle. The effective stress tensor of 
imperfect fluid under cosmological principle now reads 

(Tij = -Sij(p - (d k v k ) (6) 
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In the general theory of relativity, the explicit expression of moment-energy tensor then reads 

Tfiu = pU^Vu + (p - 3FC) V ( 7 ) 

where p is energy density, p is isotropic pressure of perfect fluid, U 11 = (1,0,0,0) is the four-velocity of the cosmic 
fluid in co-moving coordinates, h^ v = + U^U V represents the projection tensor, and 3H = d^U^. The similarity 
could be found between the classically non-relativistic expression ([5]) and the generally relativistic description about 
the applications of bulk viscosity in dynamical equations as we presently know. 

The characteristic explanation of viscosity could be classified into two types: the dynamical viscosity (bulk and 
shear viscosity) and the kinematic viscosity which is more usually called turbulence, and is normally represented by 
a characteristic Reynolds number Re. It is well known that lots of observations support the existence of turbulence, 
but it is not possible to define a universal Reynolds numbers which characterizes any large-scale complex system. 
Even though some approximation has been taken in modeling that (l2j , we should realize that kinematic viscous fluid 
models could hardly be accepted without considering the back-reaction effects. Meanwhile, the dynamical viscosity is 
naturally defined to have the feature of unification, which leads to the focus of reaches in the effects of bulk viscosity 
on cosmological evolution. In this article, we specifically explore the cosmological effects caused by bulk viscosity 
which is temperature related from molecular-like collision. 

C. approaching to viscosity models 

The differences between viscous fluid cosmological models are not only due to the different types of viscosity, but 
also to the different treatments toward the composition of the universe. So far as we know the components highly 
established both theoretically and observationally of the universe are non-relativistic dust matter (or more specifically, 
the baryons and dark baryons), cold dark matter mainly, radiation (including relativistic neutrinos and CMB) and 
the mysterious dark energy. In imperfect fluid models discussed by previous literature researches, either some parts 
or all of the components of cosmic energy-matter compositions in the universe is assumed having viscous feature. 
In this present work we take the assumption that viscosity should be a global feature which is more natural and 
rational as compatible with the cosmological principle, and try to find out whether it is possible to explain dark 
energy phenomena, which is commonly described by an exotic scalar field such as the quintessence or k-essence, or 
maybe partly caused by an observational effect of the cosmic media viscosity or due to the gravity modifications to 
the general theory of relativity at large scales. 

According to Eq. ([7]), the effective pressure which includes the contribution from bulk viscosity reads 

p = wp- 3H( (8) 

where £ is the bulk viscosity parameter, and w is the EOS parameter of cosmic fluid except for the viscosity. For 
unified fluid model (we neglect radiation since it contributes little to the late-stage evolution scenario), it is assumed 
that the universe only consists of pressless dust (baryons plus cold dark matter). The EOS parameter w in Eq. © 
should consequently be zero, and then the cosmological expansion is caused mainly by the negative pressure —3H(^ 
provided by viscosity only. 

Much efforts have been made to study imperfect models with simply constant viscosity parameter in previous 
studies. One drawback of that assumption is that it has overemphasized the differences between perfect fluid models 
and imperfect ones. We believe that the perfect fluid universe should be an evolutional approximation of imperfect 
fluid universe when it goes infinities, and the evolutional behavior of viscosity should be functioning. 

In this article we try to extend the concept which is one possible origin of cosmic viscosity — molecular-like (galaxies 
and galaxies collision, for example) interaction — to cosmology researches. On large scales we assume, that most of 
the visible cosmic matter is effectively in a thermal equilibrium state characterized by a certain temperature, and the 
interaction between galaxies or matter clumps could be approximately described in a statistical way. 

For imperfect fluid, whose viscosity is generated from molecular-like interaction and characterized by macroscopic 
temperature of the environment, the bulk viscosity can be expressed by a function £(T) as we have just discussed above. 
The CMB spectrum indicates that the universe has reached a thermal equilibrium state with a typical temperature, 
which means that the collision or interaction in the equilibrium state can be represented by the average temperature 
of ambient photons which pass through the universe and frequently interact with all luminous matter. 

Under those considerations, we can put forward in another way that global viscosity is naturally introduced by the 
gross grain approximation and related to the CMB temperature at various stages (which will be discussed briefly in 
the next subsection). 

In this paper, we adopt two appropriate molecular-like collision/ interaction models for cosmology study, from 
which the viscosity is physically generated; one is based on the Chapman-Enskog formula [32j for dilute fluid which 
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reads 

C = 2.67xl0- 6 4^ (9) 

where a is collision diameter, uj is collision integral, T* = kT/e with the maximum energy of attraction e divided by 
Boltzmann constant k, and m is molecular like mass. 

The other one is the Sutherland's formula [HJ which reads 

C = C '(^L)I^±^ (io) 

where S is Sutherland's constant, the prime represents reference values here. 

The above two models are both semi-theoretical or empirical and with first order approximations. The Chapman- 
Enskog relation assumes that molecules possess only translational kinetic energy, while Sutherland's formula assumes 
molecules are smooth rigid elastic spheres surrounded by fields of attractive force. Both forms are available for building 
cosmological models under proper simplification, since the empirically determined parameters are impossible to be 
given precisely on cosmic large scales. At low temperature (T <C 300K), a general equation which can approximately 
represent the main feature of Chapman's and Sutherland's formulae reads 

C = CoT Q (11) 

where Co is considered as a simplified coefficient, and a equals \ or | which are chosen specifically in order to 
physically represent Chapman's or Sutherland's limit formula respectively. Theoretically, a could not reach | exactly 
in Sutherland's formula (|10[) . so we pick this value as the upper limit. Thus, the Chapman's model and Sutherland's 
model of viscosity origin can be unified by Eq. (jTTJ) and extended at cosmic large scale also for late-stage evolutional 
era, while variations are provided mainly by the parameter a. 



D. The red-shift related CMB temperature 



The telling-tale CMB temperature, which is feasible to be measured at high red-shift in distant molecular clouds 
via absorption spectra Q , has provided us with lots of useful information on the universe evolution history. According 
to the current theories and observations, we are confident in assuming the linear relation between CMB temperature 
and its corresponding red-shift [(| holds: 

T CMB (z)^T Q (l + z) (12) 

where To « 2.73 the ambient CMB temperature measured at present. As we have discussed above, it is reasonable 
to combine Eq. (fT2"j) with the viscosity formulae © fTUj) in order to extend those concepts of viscosity from generally 
classical physics to cosmological physics. 



II. MODELING AND CONSTRAINING 



A. unified fluid model 



1. basic equations 



According to the discussion above, the general formula adopted as our temperature related viscosity model reads 

C = CoT a (13) 

where Co is a simplified parameter and parameter a can be chosen as either i (for Chapman's relation case) or | (for 
Sutherland's formula limit case). In addition to the expression of viscosity, unified fluid model and linear relation 
between CMB temperature and the red-shift are considered together in this subsection. Now we can rewrite the 
mathematical expressions of those two premises as: 



P = -3#C 
T(z)=T (l + z) 



(14) 
(15) 
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By combining the two Einstein's FRW equations, i.e: Eqs. ([T] [2]) with our temperature related viscosity model: 
Eqs. (fI5 ]) -([T5 1) we have 



-H 2 + H = 127rGH( T*(l + zf 



(16) 



It is easy to convert the differential equation with respect to cosmic time into that with respect to red-shift z or scalar 
factor a, by utilizing 



The general solution reads 



E(z) 



e -f(i+*r 



(17) 



(18) 



where (3 = 127tC ^° t o w jth Hq represents the Hubble parameter at present. In order to nail down the integral constant 



C, we set the boundary or initial condition E(z = 0) = 1 , so that C 



-Ei[—£-]. The symbol Ei represents 



the exponential-integral-function Ei(x) = — ^r-dt. We briefly show the behavior curve of the Ei(x) function in 
Fig. ©. 




FIG. 1: The curve of function Ei(x). A turning point occurs when x = 0, and the left branch (Ei(x < 0)) converges to zero as 
x goes to minus infinity. Since — — (1 + z) a is always negative if we ask the redshift not less than -1, the value of Ei in Eq. (|18[) 
is confined under zero mainly. 

Since the lower limit of viscosity is zero (viscosity free) or a constant (cosmological constant-like case), parameter 
P should be positive and so is the value of — j-Ei[— @- (1 + z) a ] when we consider the recent evolution of the universe, 
say from z=0.4 to the future. It is possible that Eq. (|18[) may suggest a mathematical and physical turning detour 
at the point C = ^Ei[— ^ (1 + z) a ], and this case only happens when the value of C is negative (also see Fig. JT])) 
according to the discussion above. So the negative C may denote bouncing universe models while positive C indicates 
models without the bouncing solution character. If C equals to zero, a will be related to /3 by the boundary condition, 
and consequently with a strict control on the freedom of parameter j3 to make the model less arbitrary. 



2. observational data constraints 



To confront the models for extreme cases with the observational data-sets we use, for convenience, the reduced 
dimcnsionlcss Hubble parameter h = Hq/(100 kms~ 1 Mpc~ 1 ) at present with the parameter j3. Since the matter- 
energy densities of dust matter and dark energy are coupled together in the unified fluid model, the observational 
H(z) and SNc la data sets (will be briefly discussed in the appendix) could be expected to give appropriate constraints 
on our model. 

When the model with case a = i, the observational H(z) and SNc la data sets together generate an overlap region 
with la and 2a confidence ranges for the parameter pair (/3, h) ( results are shown in the left panel of Figs. @). 
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FIG. 2: The 1 and 2 a confidence ranges of (ft,h), constrained by the observational H(z) and SNe la data sets. The left panel 
shows an ideal overlapped area for model with the a = 1/2 case, while to constraints for the extreme model case a = 3/2 the 
data-sets fittings from observational H(z) do not overlap with those from the union 2.1 of SNe la data-sets. 



While at case a = §, the H(z) and SNe la data-sets do not give overlap region for the parameter pair ((3,h), ( 
results are shown in the right one of Figs. ©). Recall that when we first mentioned the case with a = |, we take it 
as an extreme limit of the model parameter a in Sutherland's formula. Now we see that the limit could not only be 
reached theoretically, but also unfavored by observational data fittings. 

We find out that the unified fluid models are with interesting features as the dark energy phenomena required. When 
the parameter a < 3/2, the unified fluid model and the theoretical explanation of viscosity origin (the molecular-like 
interaction theory) are acceptable by the data fittings as shown. We keep in mind that so far by the global fitting any 
feasible cosmology model to include the dominated dark energy should not deviate the ACDM model away obviously. 
That conclusion leads us to move on to the study of effective viscosity functioning on the basis of the well-known 
ACDM model. 



B. ACDM-V model 



1. basic equations 

In this subsection, we focus on the application of viscosity to the ACDM model. The expression of viscosity still 
follows Eq. (fTTj) which contains both the physical content of Chapman's relation and Sutherland's formula, while the 
corresponding Einstein's component equations which include the constant dark energy A term in the flat FRW metric 
read 

^= 8 -f P = 8 -^( Pm + PA ) (19) 

H 2 + H=-^(p + 3^) (20) 
3 c z 

where p m is energy density of pressureless matter, pa represents energy density of the constant dark energy, and c is 
the light speed in the vacuum (in natural units c equals to one). The expression of the effective pressure p including 
the viscosity contribution in Eq. (|20p then reads 

P = -Pac 2 - 3H( (21) 



where the first term is the effect of constant dark energy which EOS parameter wde = — 1 exactly, and the second 
term comes from the effect of viscosity. We do not regard viscosity as one peculiar component of cosmic matter-energy 
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in the universe, but as a gross grain approximation effect of the thermal and dynamical effect of interactions in the 
fluid which permeates the universe. 

The combined form of Eqs. (|T9| (|2T|) reads 



which can be parameterized as 



H = --fi m0 ff (l + *) 



12ttGH( 



(22) 



dE Al+z 

— A 2 

dz E 



-— (1 + z)^ 



(23) 



where A — §f2 m o and B = 12ttG(qTq with a equals to 



5 or I 



according to the specific formulas of viscosity. The 



power index ±i in the above equation corresponds to Sutherland's formula and Chapman's relation respectively. 
We name this model ACDM-V where V is short for the viscosity. For convenience, letter "a" and "b" are added as 
postfixes when a equals to ^ and | respectively. The new model differs from the standard ACDM with an additional 
term -^-(l + z)*^ in Eq. (|23|) . which has no clearly analytical solution. We will use numerical method by data statistic 
analysis in model constraining. Notice here exist two theoretical restrictions for Eq. (f2"3"]l . i.e: E{z = 0) = 1 which is 
a natural boundary condition, and B > which must be ensured due to the physical existence of viscosity. 



2. observational data constraints 



In this subsection we have proposed a model which explicitly announces that the universe consists of baryons, cold 
dark matter, and constant dark energy with viscosity, and both the dark energy and viscosity contributes to the 
effective pressure. In addition to the observational H(z) and SNe la data sets, the baryon acoustic oscillation (BAO) 
data set which has strong constraint on dark energy models will also be included in model constraining. (Applications 
of BAO data will be briefly discussed in the appendix and we should mention that several observation H(z) data 
points in the general H(z) data catalogue are from BAO. We will not treat them twice in doing the statistic analysis 
[18j for the models.) 

The best-fit results for model Va and Vb are listed in Table (jj), besides which the confidence ranges of parameter 
pair (r2 m o, B) of both models are shown in Figs. ([3]). 



TABLE I: Best-fit parameters of ACDM-V models by various data sets 

data sets applied H(z) SNe H(z)+BAO 
model- Va 

xLn 15.058 562.225 16.621 

Q m0 0.367 0.283 0.284 





= UttGCot}) 31.422 


0.870 


4.533 


Ho 


{km/s/Mpc) 73.332 


69.884 


70.549 


model- Vb 




Xfnin 15.037 


562.224 


16.740 




n m0 0.589 


0.290 


0.288 


B( 


= 127rGCoT!) 58.902 


1.610 


1.921 


H 


(km/s/Mpc) 73.748 


69.885 


69.882 



The observational H(z) or SNe la data alone does not provide much tight constraints on parameter pair (Q m o,B), 
until BAO data is taken into consideration. The overlapped 2a confidence ranges in Figs. §3§ favour positive value 
of B, and confine it under the order of the magnitude of 10, while the overlapped confidence ranges of £l rn o are also 
acceptable, which lie between 2.5 and 3.0. The confidence ranges of H are not shown here, since the constraints 
on the Hubble parameter by each individual data set are tight and give out similar values around 70 km/s/Mpc. 
Besides, among the three parameters, the values of B and fi m o matter most to our newly viscous cosmology model 
constraining. 
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FIG. 3: The la and 2a confidence ranges for parameter pair (Q m o,B), constrained by the observational H(z) (dashed line), 
SNe la (solid line) and BAO+H(z) data-sets (thick solid line). The contours of confidence ranges given by the observational 
H(z) and SNe la data-sets exceed the figures' boundary, but we are mainly interested in the overlapped area. Constraints from 
the observational H(z) data alone deviate a little from those of the other two, which indicates that the observational H(z) data 
lacks enough precision to be an independent data set (unlike the good quality of SNe la data-sets as the Union 2.1) for model 
constraining. But considering the future astrophysics observations progress the observational H(z) data is still our first choice 
when a quick fitting or an auxiliary data set is required. 



III. DISCUSSIONS ON THE ACDM-V MODEL 



In the SI units, the expression of viscosity reads clearly as 

where a equals 1/2 in model- Va and 3/2 for extreme limit in model- Vb, and c = 2.99 x 10 8 m/s is the light speed in 
vacuum. The above equation can be expressed explicitly as 

C« 1.151 x B(l + z) a x 10 6 Pa- s (25) 

The evolutional behavior of viscosity is related to the expansion scalar a = j-j-j which has been pointed out in 
Eq. (|25|) . The magnitude of viscosity will be gradually reduced by the on-going expansion. This provide us an 
physically acceptable transition mechanism of fluid from imperfectness to perfectness during the cosmic evolution. 

According to the results given by model constraining, both model Va and Vb are acceptable. The results provide 
an allowed range for the value of £?, and consequently set an upper limit for viscosity parameter. Since in this present 
study we are mainly interested in the evolutional effect of viscosity on late stage cosmological evolution of the universe, 
in the following we will show how much deviation the viscosity model constrained by the background data-sets may 
provide against the standard scenario, the ACDM model. 

For convenience, we re-express the ACDM model in a simple form with the f2 m o as the cosmic observational matter 
composite fraction with repeat to the critical density of the observable universe today 

HIcdm = H$[n m0 (l + zf + (1 - O m0 )] (26) 

According to the latest WMAP-9yr observational results we set the parameters f2 m o = 0.287 and Ho = 69.32 by 
the global best-fit results for the standard cosmology model. In order to observe the effects of viscosity, we also set 
the parameter A = 0.427 and Hq = 70 for model Va and Vb, while the value of parameter B varies from to 10 
roughly according to the previous fitting results. 

The contrasts with H-z relation between the viscous models and standard cosmology models are illustrated in 
Figs. @. Viscosity contribution suppresses the value of Hubble parameter at high red-shift, which means that it 
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increases the cosmic age. As we know today that the age of the universe is 13.772 ± 0.059 Gyr as predicted by the 
standard cosmology model with observations [l[. For the possible largest value of parameter B, the standard cosmic 
age could be extended by 1.012 and 1.041 times by Va and Vb models respectively. Generally speaking, the age of 
the universe predicted by the ACDM-V models will not be older than an upper limit as 15 Gyr. 

We are also interested in the fate of the universe predicted by the viscous models. To some extent, the ACDM-V 
model can be regard as a phantom- like dark energy model according to Eq. (j2"Tj) . (also see Figs. ([5])). As the expansion 
goes on the observational universe will tend to infinities, the cosmic media viscosity will fade out or be negligent totally 
and the imperfect fluid will gradually become perfect one to satisfying accuracy. 

There are basically four types of future singularity [tJ. Type I (Big Rip): the scale factor becomes infinite at a 
finite time in the evolution future. Type II (Sudden Rip): the scale factor and energy remain finite at the rip time 
when pressure becomes infinite. Type III: the scale factor remains finite when energy and pressure become infinite in 
the finite future. Type IV: the scale factor remains finite in the finite future when energy and pressure vanish, and 
the higher derivatives of H diverge. We may extend the classfication to include the EoS singularity as well. 

When the red-shift approaches —1, the Hubble parameter remains finite as shown in Figs, (jlj, which means the 
Big Rip is avoided by model Va and Vb respectively; according to Eqs. (|21[) . (fT3|) and (|T5j) . pressure will always be 
finite, so Type II and III singularities are also absent. The last kind of singularity is apparently not predicted, so the 
ACDM-V model will not experience the above four types of singularity as well as the EoS singularity obviously. 

Taking into account of the no-rip, little rip and pseudo-rip [8[ cases, there are seven types of fate for the universe 
evolution. The ACDM model certainly predicts the no-rip future. We have already known that during the expansion 
of the universe, our viscous model will gradually evolve into the standard perfect one. With the Hubble parameter 
predicted to be constant (dominated by a cosmology constant like term) when time goes to infinity, the ACDM-V 
models are free of the kind of rip destiny. 




-1.0 -0.5 0.0 0.5 1.0 1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 

Redshift z Redshift z 



FIG. 4: The H-z relation diagram with H(z) data as reference background. The thick line represents the H(z) behavior from 
the ACDM model, while the shadowed area represents the possible H(z) values of the specific model with the value of B varying 
from range to 10 for model comparings. 



IV. CONCLUSIONS 

In the above sections we have explored the possible effects of the existence of temperature related viscosity on 
cosmological evolution by proper modeling which is supported by observed data sets fittings. Similar research works 
in literature have proposed a variety of kinds of viscosity forms, from constant viscosity parameter Q fill ] to com- 
plicated turbulence phenomena related forms as possibly appearing in the complex cosmos media physics fl2l ]. and 
the temperature related viscosity in this present work is proved to be one possibility of them with suitable models 
supported by observational data sets well confronted. 

In this article, at first we have attempted to conceive the unified dark energy (matter, dark matter and dark energy 
unified) models, expected that the temperature related viscosity may be one of the better candidates for dark energy 
phenomena. The observed data-sets fitting results have shown some inconsistence between the data-sets fitting from 
the SNe la confidence ranges against those from the observational H(z) data-sets for the viscous dark energy model 
with the extreme case of a = 3/2. The not so precise observational H(z) data-sets fitting results (for this point we can 
see clearly also from the listed table one for the H(z) fitting results accordingly) for the specific model with parameter 
a = 3/2 case (actually in the Sutherland's relation, the parameter a < 3/2) imply we need further proper dark energy 
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FIG. 5: The w-w' phase space of the ACDM-V models, where the prime represents instead of the with N=lna(t) and 
w here is the effective EOS parameter since we hope to give a clear behavior of the effective equation of state especially when 
w = — 1 rather than to illustrate the dark energy model classification by Caldwell and Linder parametrizasion (Refs. The 
thick line represents the w-w' curve of the ACDM model, while the dashed line shows the limit of viscous models in phantom 
area with value B = 10. 



model building, which means this kind of viscosity dark energy could only be one possible component describing some 
properties for currently cosmic dark energy phenomena, while the model with a = 1/2 case is satisfying. Similar 
researches have pointed to basically the same results, which indicate that unified fluid models with viscosity can mimic 
the evolutional behavior, but usually cause damping of density perturbations @ (lOj . 

Following the extension of our first trial we have then proposed the ACDM-V model, which is a modification of the 
standard ACDM model accordingly, turning the idealized perfect fluid cosmology model into the practically imperfect 
one with proper viscous energy-momentum tensor functioning and is fitted satisfactorily with the currently main 
observational data-sets. The upper limit on parameter B given by data-sets constraints is still loose, since when B 
equals to 10 the DE density may become negative at high red-shift, but notice that the null energy condition (NED) is 
not violated. That is mainly caused by the data-sets whose observational range are limited under red-shift z = 2. Since 
late-stage evolution was mainly discussed, CMB-shift parameter has not been taken into account, which can provide a 
richer content about the past history of the universe. For consistency, we strongly suggest that a theoretically proper 
upper limit on B may lie around the magnitude of order of 0.1 in order to avoid that crossing phenomenon of DE 
energy density. Introducing interaction between DE and DM into this model may also expected to raise the upper 
limit. For illustration, the evolution behavior of DE density of ACDM-V model with B = 0.1 arc shown in Fig. ((6]). 




Vb 



-1 012 3 45 -1 012345 

red-shift z red-shift z 



FIG. 6: The density evolution of DE, when B = 0.1 as an example, without singularity or crossing behavior at least after the 
occurrence of decoupling (z < 1100). 
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In the Figs. (|5|) for numerical comparisons, we can see the parameterized equation of state (EoS) parameter evolution 
lines of the ACDM-V and the standard ACDM model can match at present stage and in the evolution future, to form a 
consistent picture which is in accordance with the no-rip future predicted by the new model. In addition to increasing 
the cosmic age and being free of future singularity, we would like to say the smooth transition mechanism is another 
pleasing feature of the ACDM-V model. In previous researches conducted in some cosmological dark energy model 
buildings, imperfect fluid models remain to be with imperfect fluid all along the evolution of the universe. In this 
present work, the model Va and Vb could converge to the cosmic perfect fluid model, the globally fitted well ACDM 
model, as the cosmic expansion goes on to infinity to dilute the viscosity effects, by providing a smooth transition 
mechanism from practically imperfect fluid models to perfect fluid ones. 

Considering about the complexities and difficulties for the physical viscous media descriptions in observational 
universe, in the present work we have not included the density perturbation discussions which may give much tighter 
constraints on model buildings, like via the integrated Sachs- Wolfe effect or matter power spectrum simulations 
(see Refs. @ [ll| QJ] [3). We strongly believe that more precise observed data sets given by future astrophysics 
observation instrumental achievements can certain provide higher capacity performance in model constraining and 
constructing and will finally pin down the physics for long puzzling cosmic dark matter problems and mysterious dark 
energy phenomena. 

The currently possible origins and explanations of viscosity physics may remind us of that trying to express viscosity 
effect in only one form could be too simple. The rich viscosity effects in cosmos may be an assemblage of a variety of 
dynamical and kinematics effects of locally cosmic motion and globally cosmos large-scale evolution, which may show 
different features when observed on different scales and /or in different evolutional stages. We may say this present 
article research reveals only a tip of the whole iceberg. Nevertheless, continuous study endeavors to the mysterious 
dark sector physics will surely shed light towards the fundamental understandings of our universe. 
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Appendix 

For the possible benefits to the concerned academic community with our hope we write out in details the data 
fitting preparations in this work. We employ heavily in this presentation the \ 2 statistical method [l6| E3 f° r 
model buildings with probability distribution function which reads (for measuring confidence range of parameters or 
coefficients) 

P{$)=Ne-^ (27) 
where TV is the normalization coefficient to ensure J P(d)M = 1. 



A. Observed Hubble Parameter data (OHD) sets (the observational H(z) data sets) 

The H(z) data-sets contains 28 measurement points so far. It provides a clear and direct picture of H(z) evolution 
over the range 0.1 < z < 1.5. OHD depicts the relation between the Hubble parameter H(z) and the according 
red-shift z, which is also related and described directly when conceiving cosmic models. The value catalogue of this 
data-set are reported in Table ((TTJ) . The x 2 calculation for OHD reads 

XOHD - />t 1 

* * /T„- 
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TABLE II: Observed Hubble parameter Data, Refs. [THIH 



z H (km/s/Mpc) error 



0.090 


69 


12 


0.170 


83 


8 


0.270 


77 


U 


0.400 


95 


17 


0.900 


117 


23 


1.300 


168 


17 


1.430 


177 


IS 


1.530 


140 


11 


1.750 


202 


40 


0.480 


97 


62 


0.880 


90 


40 


0.179 


75 


1 


0.199 


75 


5 


0.352 


83 


14 


0.593 


104 


13 


0.680 


92 


8 


0.781 


105 


12 


0.875 


125 


17 


1.037 


154 


20 


0.24 


79.69 


3.32 


0.43 


86.45 


3.27 


0.07 


69.0 


19.6 


0.12 


68.6 


26.2 


0.20 


72.9 


29.6 


0.28 


88.8 


36.6 


0.44 


82.6 


7.8 


0.60 


87.9 


6.1 


0.73 


97.3 


7.0 



B. Type la Supernova (SNe) data set (Union 2.1) 

SNe la data is believed to provide a solid piece of evidence for cosmic acceleration at present era of evolution of the 
Universe. Type la supernova are regarded as standard candles for distance measurement in observation universe, with 
same absolute magnitude AI ~ —19.3 ±0.3. Union2.1 data set which is used in this paper contains 580 measurement 
points. The theoretical distance modulus is defined as 

Lkh = m t h - M = 5log w D L (z) + fi (29) 

where Dl = (1 + z) e(z') ' un der the assumption that our universe is spatially flat, at least within the space-time we 
are able to have observed, = 0; and fiQ — 42.38 — 5logioh,, where h is the dimcnsionless reduced Hubble parameter 
at present. The x 2 fitting formula for SNe la data then reads: 

580 / \ , a\ 

^^Y f^-^^ f (30) 

Since some of our viscous models arc very complicated, we may reformulate the SNe la data set by converting the 
(i ~ z relation into Dl ~ z relation through the following equations: 



(31) 



13 



besides, we should convert the old error bar value of each measurement point into the corresponding new error value 
with 

CTnewH) = D L ( Zi ) [10^-1]. (32) 

Under this transformation of data sets, we are free from the Logio calculation while experimenting \ 2 fitting. Thus 
the reformed \ 2 calculation is 

2 _^rD Lobs {zi) - D Lth (zi]fi) 2 

~ 771 J ( ' 

. . u new \ L J 

l=l 



C. Baryon Acoustic Oscillations (BAO) data set 



The BAO is originated from perturbated sound waves in the early universe when primordial sound wave could 
travel at relativistic speed through space-time. We have gathered data from 6dFGS (2f| (Beutler et al. 2011), Wiggle 
Z [1| (Blake et a£2011), and SDSS \2J^ (Eisenstein et aZ.2005). The data sets will be listed in Tabic dm} below. The 
statistically % 2 analysis for BAO data-sets reads 

2 _ M^^ Aobsjz,) - A th {zj;d) ^ 2 
Xbao — 7_\ J > 



where 



-°^ 1/3[i f " wf" (35) 



z 



TABLE III: Baryon Acoustic Oscillation Data, Refs. [2514271] 



z 


A 


error 


0.106 


0.526 


0.028 


0.35 


0.469 


0.017 


0.6 


0.452 


0.018 
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